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Abstract 

We perform a multilepton channel analysis in the context of the Large Hadron 
Collider (LHC) for Wilkinson Microwave Anisotropy probe (WMAP) compatible points 
in a model with non- universal scalar masses, which admits a Higgs funnel region of 
supersymmetry dark matter even for a small tan f3. In addition to two and three- 
lepton final states, four-lepton events, too, are shown to be useful for this purpose. We 
also compare the collider signatures in similar channels for WMAP compatible points 
in the minimal supergravity (mSUGRA) framework with similar gluino masses. Some 
definite features of such non-universal scenario emerge from the analysis. 
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1 Introduction 



Low energy Supersymmetry (SUSY) [1] is a strong candidate for physics beyond the Standard 
Model. A general framework like the Minimal Supersymmetric Standard Model (MSSM) 
[1,2], however, suffers from a large number of a priori unrelated parameters. This lack 
of predictability can be minimized only if one assumes a definite mechanism for breaking 
SUSY. The minimal supergravity (mSUGRA) [3] model that assumes gravity to be the 
mediator between the hidden sector wherein supersymmetry breaks and the observable sector 
(where the MSSM exists) is relatively economical in this respect. Starting from only a few 
parameters, renormalization group evolutions (RGE) and the use of radiative electroweak 
symmetry breaking (REWSB) generate all the MSSM parameters at the electroweak scale. 

The mSUGRA scenario has a remarkable simplicity of principle, an economy of param- 
eters and features that at least partially ameliorate potentially disastrous consequences in 
low energy physics. From a more agnostic standpoint, however, there is no strong reason 
to restrict ourselves to such universal models. For one, even with gravity conveying super- 
symmetry breaking, the soft SUSY-breaking terms need not be universal at the supergravity 
scale, but would depend on the structure of the Kahler potential. Similarly, large non- 
universal corrections may accrue to the soft parameters as a result of the evolution between 
the Planck scale and the gauge-coupling unification scale (Mq — 2 x 10 16 GeV) [4]. These 
and other related issues have led to several studies of non-universal scalar [5-15] and gaugino 
mass [16,17] models. Non-universal scalar masses may appear due to a non-flat Kahler met- 
ric [18], or, for example, from 50(10) D-terms [12-14]. However, any such nonuniversality, 
at the electroweak scale, would lead to low-energy flavor changing neutral current (FCNC) 
processes (through SUSY loops) [19]. The existing data on flavor physics thus impose severe 
constraints on any nonuniversality in scalar masses, in particular for the first two families. 
The restrictions on the third generation scalars (and the Higgses) from FCNC data are not 
too severe though. 

It turns out that both FCNC and CP-violation constraints may be best tackled if one 
assumes the first two generations of scalars to be multi-TeV and (quasi-) degenerate in 
masses [20] 1 . Clearly, allowing universal scalar masses at the gauge coupling unification 
scale would not satisfy the above objective because either (i) the REWSB constraint would 
prohibit such large scalar masses for a reasonable set of values of the gluino masses, or (ii) 

1 We remind the reader that satisfying constraints imposed by electric dipole moments of electron and 
neutron would require very large scalar masses if we like to have finite values for the CP-violating SUSY 
phases. 
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one must have very large gaugino masses, so as to allow very large scalar masses, thereby 
worsening the fine tuning problem [21]. We recall that, within the MSSM, the naturalness 
problem and its solution revolve around the third family, as well as the gaugino and Higgs 
scalar mass parameters. As long as the third generation scalars and the electroweak gaug- 
inos are on the lighter side, any quantitative measure of naturalness would stay within an 
acceptable domain. Furthermore, constraints from FCNC and CP-violation are relatively 
weak in such a scenario with an inverted mass hierarchy [22,25]. 

In this work, we consider a particular non- universal scalar mass scenario (NUSM), namely 
that of Ref. [5]. The model addresses the FCNC issue by invoking very large masses for the 
first two generations of squarks and sleptons. As is well-known, such a solution is difficult 
to achieve within the mSUGRA scenario as the requirement of REWSB prevents the scalar 
masses from being too large. In the present context, this is circumvented by allowing the 
third generation squark masses and the Higgs scalar mass parameters to be small. This very 
smallness also serves to keep the degree of fine-tuning within control. 

As far as the third generation sleptons are concerned, a very small SUSY-breaking mass 
at the GUT scale is not phenomenologically viable since the larger Yukawa coupling serves 
to drive down the mass of the lighter stau, thereby rendering it the lightest of the super- 
symmetric partners (LSP) at the electroweak scale. Consequently, the SUSY-breaking mass 
in this sector has to be sizable 2 . Rather than introducing a new parameter, we shall assume 
it to be same as that of first two generations of squarks or sleptons. To summarize, at the 
GUT scale, all sfermion masses are diagonal; and, apart from those pertaining to the stop 
and the sbottom, are universal. The last-mentioned, along with the Higgs scalars, have a 
vanishing mass at this scale. While this construction might seem artificial, note that this 
accords a special status to only those fields that are expected to play a direct role in EWSB. 
Interestingly, the model satisfies the WMAP constraint [26] on neutralino relic density for 
a large region of the parameter space without requiring any delicate mixing of Binos and 
Higgsinos. For simplicity, we confine ourselves to a universal gaugino mass and a vanishing 
trilinear soft-breaking parameter (A ) at Mq. 

We investigate how such a scenario can leave its fingerprint on numbers measured at 
the Large Hadron Collider (LHC). Such fingerprints are of value if ways can be devised 
to distinguish this scenario from an mSUGRA one with, say, similar gluino masses. For 

2 However, in analyses with Higgs-exempt no-scale SUSY model [23] or in a model with gaugino mediation 
[24] one may avoid such charged LSPs at the electro-weak scale by using non-zero Higgs scalar masses at 
the unification scale. In these scenarios the no scale boundary conditions are also valid for sleptons. 
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this, one has to perform a multichannel analysis [15,17,27,28] studying several final states 
simultaneously. 

A promising signal of supersymmetry (with a conserved R-parity) comprises large missing 
transverse energy, accompanied by jets and leptons with varying multiplicities. An analysis 
in different channels, compared with that of a similar mSUGRA scenario may lead to a 
significant hint of the non-universality. In the present analysis, we assess the accessibility 
of our non-universal scalar mass model (NUSM) at the LHC We find that the direct pair 
production of stops and sbottoms as well as their cascading down from gluino decays lead to 
the possibility of four-lepton final states as a distinct signature of this scenario. Additionally, 
we also analyse the two-lepton and the three-lepton final states. This includes opposite sign 
dilepton, same-sign dilepton and trilepton final states. All these analyses are done also for 
mSUGRA so that the multipronged approach of analysing for different channels may become 
more conclusive. 

The paper is organized as follows. In Section 2 we describe the NUSM model, apply 
cosmological constraints on neutralino dark matter and use low energy constraints such 
as those from b — >■ s + 7 or B s — > fi + fi~. We also identify benchmark points for our 
analyses of collider signals at the LHC. In Section 3, we pinpoint our strategies for collider 
simulations and report the numerical results. Finally, in Section 4, we summarize our results 
and conclude. 

2 The Non-Universal Scalar Mass model (NUSM) and 
benchmark points 

2.1 The NUSM parameter space 

The NUSM model [5], at the scale Mq, is characterized by five parameters, namely, 

tan/3, mx/2, mo, A and sign(fi) . (1) 

The parameters, here, play roles similar to those in mSUGRA except for a subtle and im- 
portant difference in the scalar sector. Masses of the first two generations of scalars (squarks 
and sleptons) and the third generation of sleptons are assigned the value m . However, the 
Higgs scalars and the third family of squarks have vanishing mass values at Mq. Here, m 
is allowed to be up to tens of TeVs. As has already been stated, we limit ourselves to a 
vanishing Aq in this analysis. We have considered fi > in this analysis. 
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The NUSM admits a smaller pseudoscalar Higgs boson mass wia on account of the Large 
Slepton Mass (LSM) renormalization group effect [5] for large m . With such a m , the LSM 
effect causes m 2 H to become large and negative and this may happen for even a small tan/9. 
This, in turn, reduces the masses of the pseudoscalar Higgs boson (A), the CP-even heavy 
Higgs boson (H) and the charged Higgs bosons (H ). In this scenario, fi is quite insensitive 
to a change in m [5], since the Higgs and the third-generation squark masses at Mq are 
free of the latter. In fact, /i is completely independent of m up to one-loop, whereas the 
two- loop contributions to its RGEs result in only a tiny dependence on m . Recall that, 
in mSUGRA on the contrary, decreases significantly with an increase in m Q . Whereas 
this led to a very small for a large mo in mSUGRA, giving the so called Hyperbolic 
Branch/Focus point (HB/FP) [29,30] region that is close to the upper limit of satisfying 
REWSB for a given mi/2, there is no HB/FP type of effect in NUSM and \i stays reasonably 
independent of mo- It turns out that the lightest neutralino is highly Bino-dominated (with a 
small Higgsino admixture) throughout virtually the entire parameter space of NUSM. Along 
with the resonance condition, namely, 2m^o ~ (m#), the small Higgsino content allows 
the LSP to have the right degree of pair-annihilation via s-channel Higgs-exchanges, so as 
to satisfy the WMAP limits on the neutralino relic density. It is important to note that, 
excepting for LSP-stau coannihilation, the Higgs-pole annihilation mechanism is the only 
one in NUSM that reduces the relic density from overabundance to an acceptable degree 
of abundance. Thus, unlike in models such as the mSUGRA, here one does not need any 
delicate mixing between a Bino and Higgsinos in order to satisfy the WMAP data. Such 
Higgs-pole annihilations that occur for large tan/3 in mSUGRA is typically known as the 
funnel region [31,32]. NUSM has an extended funnel region that spans from low to high 
tan j3. 

2.2 Cosmological and low energy constraints in NUSM 

Assuming that dark matter was generated thermally, the limits on the cosmological relic 
density from the WMAP data [26] impose severe constraints on supergravity type of models 
wherein the lightest neutralino Xi becomes the LSP for most of the parameter space [33,34]. 
We now perform an analogous analysis for the NUSM. For a given set of parameter (vide 
eq.(l)) values, the supersymmetric particle spectrum is generated using SuSpect v2.34 [35]. 
This, then, is used as an input to micrOMEGAs [36] for computing the neutralino relic density. 
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The recent WMAP data [26] stipulates that, at the 3a level, 

0.091 < VtcDAih 2 < 0.128 , 



(2) 



where VLcdm is the dark matter relic density in units of the critical density and h = 0.71 ± 
0.026 is the reduced Hubble constant (namely, in units of 100 km s _1 Mpc -1 ). 
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Figure 1: (a) WMAP allowed regions in the mi^ — m^ plane for tan f3 = 10 and Aq = with 
fi > /or NUSM are shown in red dots. Lighter Higgs boson mass limits are represented 
by solid lines. Dot-dashed line refers to b — > limit. The entire region is allowed by 
B s — > data, (b) Same as (a) except that tan = 15. The B s — > bound is shown 

as a long-dashed line. This eliminates a small strip of region below the discarded top (cyan) 
region. (c)Same as (b) except that tan/3 = 40. 

In Fig.l, we display the allowed regions in the m^-mo plane for three values of the ratio 
of the Higgs vacuum expectation values, namely tan/3 = 10, 15 and 40. The thin (cyan) sliver 
at the bottom is ruled out as, for such values of the parameters, the lighter stau becomes the 
LSP. The upper (cyan) region is rejected primarily on account of the failure in the breaking 
of the electroweak symmetry via radiative means. In other words, for such parameter values, 
m 2 A does not acquire a positive value through RG flow. Close to the boundary of this region, 
several other phenomenological constraints become important. The most important of these 
pertain to (i) the LEP2 and Tevatron lower bounds for sparticle masses, (ii) sfermions 
turning tachyonic, or (m) the appearance of charge and color breaking (CCB) minima. To 
be allowed, a parameter point must evade all these and other such constraints. Specific 
details may be found in Ref. [5]. 



6 



Highlighted (in bold — red — dots) in Fig.l are examples of parameter points that satisfy 
the WMAP data. There are two distinct regions with acceptable relic density as already 
mentioned in Sec. 2.1. (a) The Higgs pole annihilation region (also known as the funnel 
region) is characterised by 2m^o ~ VfiAyVfin- In this particular scenario, it extends over the 
full range of m^j under consideration. The Higgs pole annihilations may occur through 
s-channel pseudoscalar Higgs boson (A) or CP-even neutral H or h-bosons. NUSM has a 
bino-dominated LSP similar to what occurs in mSUGRA in its funnel region that satisfies 
the WMAP data. Similar to the case in mSUGRA, the WMAP satisfied parameter regions of 
NUSM is also dominantly characterized by the pseudoscalar Higgs boson mediated resonance 
annihilation. The exact or near-exact resonance regions have very large annihilation cross 
sections resulting in a high degree of under- abundance of dark matter. The resonance region 
that satisfies the WMAP data may be a few Ta/h away from exact resonance. The widths 
Ta/h (of A/H bosons) can be fairly large {e.g. Ta/h ~ 10-50 GeV). The WMAP satisfied 
regions fall on either side of the exact resonance condition thus showing two branches in the 
figure. (6) The second region, just above the lower ruled-out part, corresponds to the case 
where the lighter stau is nearly degenerate with the LSP, leading to very efficient LSP-stau 
coannihilation, thereby reducing the relic abundance to acceptable levels. 

Also imposed on Fig.l are the pertinent low-energy constraints. Whereas non-observance 
at LEP2 impose a strict bound of 114.4 GeV on the SM Higgs [37], with recent negative 
results from Tevatron [38] ruling out even somewhat heavier Higgses, the translation of this 
bound to the MSSM case needs careful consideration. Apart from the parameter-dependence 
of the cross-sections at LEP2/Tevatron, one needs to account for the uncertainties in com- 
puting the mass of light Higgs boson [39], originating primarily from momentum-independent 
as well as momentum-dependent two-loop corrections, higher loop corrections from the top- 
stop sector etc. Numerically, this amounts to about 3 GeV, and we have taken that into 
account in drawing the solid lines representing this constraint. Additionally, a part of the 
NUSM parameter space is associated with very light rriA even for a small tan (3 and this 
may lower the lighter Higgs boson lower bound to a value much smaller than that of the SM 
Higgs boson limit. We will revert to this while discussing the NUSM benchmark points. 

A low energy observable of particular importance is the decay rate for b — > 57 rate [40,41], 
which, at the 3cr level, reads [42] 

2.77 x 10^ 4 < Br(b s 7 ) < 4.33 x 10~ 4 . (3) 

We used micrOMEGAs [36] for computation of b — >• 57 that in turn refers to Refs. [40,41] for 
actual computation. Typically, b — > s 7 disfavours the small m^j region where the rate is 
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below the lower limit. Note, however, that the usual estimation assumes a perfect alignment 
at high energies between the quark and squark mass matrices. In other words, the (super- 
)Cabibbo-Kobayashi-Maskawa matrix operative for supersymmetric diagrams is assumed to 
be identical to the usual CKM matrix. However, if one relaxes the above assumption and 
considers even a moderate amount of b — s mixing at the GUT scale, Eq.3 is no longer an 
effective constraint for high scale models like mSUGRA. This, on the other hand, will not 
cause any significant change in the sparticle mass spectra or in the flavor conserving process 
of neutralino annihilation. We refer the reader to Refs. [43,44] for further discussions on the 
amount of model-dependence in computing Br(b — > sj) in this context. 

Since the NUSM scenario may contain a light pseudoscalar Higgs, it is necessary to 
consider the constraints from B s — > Within the MSSM, the above branching ratio is 

proportional to and tan 6 /? [45]. The recent CDF [46] limit for Br(B s — > is given 

by 

Bi{B s < 5.8 x 10~ 8 . (4) 

The branching ratio of B s — > fi + fi~ is evaluated by using micrOMEGAs [36] that in turn 
implemented Ref. [47] for the computation. The computation involves inclusion of loop 
contributions due to chargino, sneutrino, stop and Higgs exchanges. The upper limit of this 
branching ratio is shown in dashed lines in Fig.l. The white regions above the dashed lines 
in Figs. 1(6, c) are thus discarded. As mentioned in Ref. [5], the intense coupling region of 
Higgs bosons that appears when rriA is very small is also ruled out in NUSM for the same 
reason. 

2.3 A few benchmark points 

The NUSM has a large volume of allowed parameter space, especially because REWSB does 
not prohibit rrio from assuming a very large value. We focus here on a few characteristic 
parameter points that satisfy WMAP as well as low energy constraints. As seen in Figs.l(a- 
c), for a given m^, the upper limit on m decreases with an increase in tan/3. As an 
example, for mi/2 = 1 TeV, mo may well be as large as 7 TeV for tan/3 = 10, 6 TeV for 
tan /3 = 15, and 1.6 TeV for tan (3 = 40. 

Here, we have preferentially explored those regions in the parameter space which give 
distinctly different low energy and cosmological signatures as compared to mSUGRA. As we 
have mentioned before, the Higgs funnel region for mSUGRA is found only for large values 
of tan /3. The NUSM is characteristically different from mSUGRA in the sense that funnel 
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Table 1: NUSM Benchmark points A, B and C (masses are in GeVs). The first five param- 
eters define the model, while the rest are predictions. 

regions exist even for small tan/3. Hence, we choose to explore one benchmark point with 
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small tan /3. Additionally, we examine a benchmark point for large tan f3. The next point to 
note is that NUSM, typically, has heavier spectra for the first two generations of scalars and 
the third generation of sleptons. Sfermions become heavier with increase in m\n. Therefore, 
if we like to probe the model in the early phase of the LHC, we would rather select m\i% to 
be relatively small for the benchmark points. 

Table 1 lists three benchmark points for NUSM. Point A (for tan/3 = 10, A = 0, my^ = 
270 GeV, m = 2.05 TeV and sign(fi) = 1), is associated with reasonably small masses for 
stop, sbottom, charginos as well as neutralinos. Furthermore, it has a light Higgs spectrum. 
All these are promising from the viewpoint of early LHC results. Point B of Table 1 refers to a 
special parameter point (tan /3 = 15, Aq = 0, vn\n = 255 GeV, mo = 2.0 TeV and sign(fi) = 
1), for which the Higgs sector is not in the decoupling [48,49] region. Thus, here we obtain 
a reduced lower limit for (close to Mz). Point C (tan/3 = 40, A = 0,m^2 = 540 GeV, 
rriQ = 1.25 TeV and sign(fi) = 1) represents a relatively heavier spectrum. However, the 
relevant parts of NUSM spectra still remain within the LHC reach. We point out that we 
have relaxed the b — > 57 constraint for points A and B. This is in keeping with the discussion 
in the paragraph following Eq.3. However, with a small displacement of the parameter point, 
we would be able to respect the constraint at the cost of having a benchmark point with an 
upwardly shifted spectrum. Both points A and B obey the constraints from B s — > 
Point C, on the other hand, satisfies all the constraints, namely, those from B s — > fi + n~ and 
b — > sj, over and above those from the WMAP data. 

We also study the collider signatures for mSUGRA scenario at points with the same 
(or very similar) gluino mass and tan/3 corresponding to the each of points A,B,C. These 
have been denoted by mSUGRA-A, mSUGRA-B and mSUGRA-C. We must mention that 
the requirements of obeying the stringent WMAP data as well as the lower bound of the 
lighter chargino mass did not allow us to choose exactly identical values of the masses of 
the gluino in each case of the mSUGRA points. This is particularly true for mSUGRA-A 
and mSUGRA-B that fall in the HB/FP zone. The high scale parameters as well as the 
low scale soft masses for these points are listed in Table 2, all of them being consistent with 
the constraint from WMAP. sign(fi) is taken to be positive and the trilinear coupling Aq 
is taken to be zero, as mentioned earlier. The corresponding low-energy spectra have also 
been generated via SuSpect v2.34 using two-loop RGEs. Full one-loop and the dominant 
two-loop corrections to the Higgs masses are incorporated. We have used the strong coupling 
a-s(Mz) = 0.1172 for this calculation, adopting the default option in SuSpect. We have 
assumed the top quark mass to be 172.7 GeV throughout the analysis, and no tachyonic 



10 



t~\ n v o wi nr nt 1 

pdldlilclcl 


mSUGRA-A 


mSUGRA-B 


ma u vji\,iT.-o 


ran p 


1 n 
1U 




1 k 

10 




/in 

4U 


m l/2 


zoo 




O^O 

zoz 






m 


07/I n 
z f 4U 




OQnn 




ZDoU 


A 


A 

U 




U 




n 
u 


sign{ii) 


1 

1 




-1 
1 




1 
1 


11 


1 on 
139 




1 
13o 




zoo 


m- g 


*7A n 
(41) 




I ZO 




1 07n 
Iz /U 




Z740 




ooon 
zozU 




osi n 
ZolU 


m h 


looo 




1 on 1 
1391 




1 7f;n 
1 ( OU 


m h 


ZZ08 




1 one 
1898 




01 on 
Z1ZU 




ZZOO 




1 on k 
1890 




01 Qn 

ZloU 


m b 2 


z / oU 




noon 

zz8z 




o/i nn 

Z4UU 


m e L 


070.1 
z7ol 




zzy4 




ZDoU 




071 /I 

z I 14 




00c 
zzOO 




ZZ 1 u 


m xf 


114 




1 1 Q 

llo 




Zoo 


m xf 


OC c 
ZOO 




zOl 




4o4 


m x% 


oc c 
zOO 




zOz 




A Q/l 
4o4 


m x% 


ICO 

loz 




1 /i n 
149 




077 
Z 1 1 




13o 




1 O A 

134 




zo / 


m x\ 


82 




O 1 

81 




1 nfi 


m A 


Z7U4 




001 

ZZ1Z 




1 7on 
1 / zU 


m H + 


2706 




2214 




1 70n 
1 ( zu 


m h 


118 




118 




119 




0.128 




0.120 




0.092 


Br{b — > S7) 


3.62 x 10" 


-4 


3.57 x 10" 


-4 


3.42 x 10" 4 




3.12 x 10" 


-9 


3.11 x 10" 


-9 


3.02 x 10~ 9 




4.60 x 10~ 


11 


1.14 x 10- 


-10 


2.56 x 10~ 10 



Table 2: mSUGRA Benchmark points A, B and C (masses are in GeVs). The first five 
parameters define the model, while the rest are predictions. 

sfermion mode has been allowed at any scale. We now comment on the differences in spectra 
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between the NUSM and mSUGRA benchmark points. For reasons that have been already 
stated, the high scale scalar mass parameters need to be chosen differently in the two cases. 
Consequently, the value of /i in NUSM is larger than that in mSUGRA, simply because the 
mSUGRA benchmark points are within or very close to the HB/FP zones. We must note 
that there is no HB/FP like effect in NUSM that would reduce fi. As a result, the chargino 
and neutralino masses in mSUGRA benchmark points are smaller than their counterparts 
in NUSM. 

Finally, in regard to the mass of gluino it is important to clarify the role of radiative 
corrections in the benchmark points of the two scenarios namely mSUGRA-z and NUSM's 
point i, where i = A,B,C. Radiative corrections comprising of gluon-gluino and quark-squark 
loops may be estimated as in Eq.5 [50]. 

m- g = m 3 (Q 2 ) + ^m 3 (5- 3 In 

01 s „ f A / 2 2 2 



q=u,..,t 

a. 



Bi(m 3 ,m q ,m qi ) +Si(m 3 ,m g ,m^ 



+ Yl 4^ m i sin ( 2 ^) Re [ 5 oK, m 2 q , m?J - B (ml, m 2 , m?J] . (5) 

q=t,b 

The Passarino-Veltman functions Bq, B\ and further useful details may be seen in Ref. [50]. 
The choice of the scale Q is not unambiguous and, in general, is defined by an appropriate 
mass scale in the theory. In SuSpect, for example, this is set equal to the geometric average 
of the values of the two stop squark masses. As we can see from Table 3, this average varies 
widely between mSUGRA-i and the corresponding NUSM's benchmark point i. As a result 
the running mass m 3 (Q 2 ) for mSUGRA-z is smaller where Q is higher compared to point % of 
NUSM, where the corresponding scale is smaller because the masses of the third generation 
of squarks in NUSM are quite smaller 3 . In general, for the given benchmark points under 
consideration, a point mSUGRA-i has a smaller running mass m 3 (Q 2 ) but has a much larger 
contribution from radiative corrections {vide Eq.5) compared to the corresponding point % 
of NUSM. We note that the radiative correction amounts that arise from gluon-gluino and 
quark-squark loops are quite different in the two scenarios. With a heavier average squark 
mass, a benchmark point mSUGRA-i has a much smaller contribution from quark-squark 
loops compared to that of NUSM point i. On the other hand, the logarithmic term in Eq.5 
is such that for mSUGRA- z the term is negative because of the fact that m 2 (Q 2 ) < Q 2 
owing to a heavier average SUSY mass scale. This leads to a large contribution from the 



i m 3 (Q 2 ) increases with a decreasing Q: see for example Ref. [51]. 
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GeV 


GeV 


GeV 


GeV 

v .l t V 


GeV 




A 


2050 


270 


357 


633 


709 


12% 


mSUGRA-A 


2740 


253 


1917 


548 


740 


35% 


B 


2000 


255 


329 


602 


674 


12% 


mSUGRA-B 


2300 


252 


1619 


551 


725 


32% 


C 


1250 


540 


903 


1197 


1278 


7% 


mSUGRA-C 


2680 


490 


1921 


1051 


1273 


21% 



Table 3: Running mass, radiative correction in percentage and pole mass of gluino in 
mSUGRA and NUSM benchmark points. The scale Q refers to geometric mean of stop 
squark mass values. 

second term of the same equation for mSUGRA-z. This, however is not true for NUSM 
where one has m\(Q 2 ) > Q 2 owing to a lighter average stop mass or a lighter SUSY mass 
scale in general. NUSM points also have significant amount of quark-squark contributions 
for the same reason. The final effect is such that smaller values of m^{Q 2 ) are overrun by 
radiative corrections in mSUGRA-A and mSUGRA-B leading to larger values of the pole 
masses nig in comparison to the values of nig for NUSM benchmark points A and B. The 
point mSUGRA-C and NUSM point C are quite competing in the above effects because of 
a larger associated and the values of nig are thus close to each other. 

3 Collider Signatures 
3.1 The general strategy 

The collider signatures, and hence the optimal search strategies, of the NUSM would natu- 
rally depend on the particular point in the parameter space that nature may have chosen. 
Rather than attempting a general, and hence non-optimal, analysis, we choose to illustrate 
the various features, concentrating largely on the three representative points identified in the 
preceding section. To start with, we summarize, in brief, the generic simulation procedure 
that has been adopted here. The spectrum generated by SuSpect v2.34 as described earlier 
is fed into the event generator Pythia 6.4.16 [52] through a standard SLHA [53] interface for 
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the simulation of pp collisions with a centre-of-mass energy of 14 TeV. 

We have used the CTEQ5L [54] parton distribution functions, the QCD renormalization 
and factorization scales both being set at the subprocess centre-of-mass energy Vs. All 
possible SUSY processes and decay chains consistent with conserved i?-parity have been kept 
open. We have kept initial and final state radiations (ISR/FSR) on. The effect of multiple 
interactions has been neglected though. We, however, take hadronization into account using 
the fragmentation functions built into Pythia. 

In Table 4, we list the total supersymmetric particle production cross-sections for each of 
the benchmark points. Also listed are the individual cross sections for some of the important 
processes, namely, gg, ti(2)^*(2) ana - ^1(2)^1(2) ana - processes with at least one chargino or 
neutralino denoted by "Xi/XiV'- We note that, for points A and B, dominant production 
accrues from stop pairs, while for point C, no production process dominates overwhelmingly. 
The other important processes include associated stop and sbottom production along with 
gluinos as well as charginos and neutralinos. That the total cross-section is much smaller for 
point C, compared to the other two, is easy to understand as the spectrum is relatively heavier 
in this case. It should be noted that while the mSUGRA and the NUSM benchmark points 
are quite similar as far as the gluino-pair production or the total supersymmetric particle 
production cross sections are concerned, they differ markedly in the dominant production 
modes. For the mSUGRA points, it is the lighter neutralinos and charginos that dominate 
overwhelmingly, whereas for the NUSM points, this role is usurped by stop-pairs and sbottom 
pairs. 



mSUGRA 


NUSM 


Point 


Total 


Xi/Xl,2 


99 


Point 


Total 


titi + ^2^2 


616; + b 2 b\ 


99 


mSUGRA-A 


11.86 


10.67 


1.18 


A 


12.42 


6.77 


1.73 


1.28 


mSUGRA-B 


12.49 


11.18 


1.25 


B 


19.92 


11.73 


2.79 


1.78 


mSUGRA-C 


0.62 


0.59 


0.02 


C 


1.23 


0.09 


0.02 


0.07 



Table 4: Total supersymmetric particle production cross-sections (inpb) as well as the leading 
contributions for each of the NUSM and mSUGRA benchmark points. 

Before we discuss the signals, it behoves us to briefly discuss the major decay modes 
(see Tables 5 and 6), for the structure of the cascades would determine the final state 
configurations. Starting with the major produce, namely the stop, for each of points A and 
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Decay modes 
(squark/gluino) 


A 


B 


C 


q — > bib 


31.0 


33.0 


28.0 


q — >• bob 


26.0 


26.0 


20.0 


g -> tit 


22.0 


21.0 


29.0 


g^t 2 t 


21.0 


20.0 


23.0 


bi — > y% 


8.0 


13.0 


7.0 


bi ->■ x% 


42.0 


52.0 


24.0 


61 -»■ xf* 


11.0 


0.0 


40.0 


61 -^tiW" 


37.0 


33.0 


3.0 


ti — >• vtfr 

-L /VI 


100.0 


100.0 


33.0 










21.0 


*i -> X°it 


0.0 


0.0 


23.0 


h -> xEt 


0.0 


0.0 


13.0 




69.0 


71.0 


0.0 


-> x?« 


10.0 


10.0 


0.0 




20.0 


18.0 


0.0 


X2 ->■ Xih 


0.0 


0.0 


90.0 


xt -+ x\w + 


100.0 


100.0 


100.0 



Table 5: The branching ratios (%) of the dominant decay modes of the gluinos, squarks and 
lighter electroweak gauginos for NUSM for the different benchmark points. 

B, it has almost a 100% decay branching fraction to b and xt- The Xi decays, in turn, 
into W ± and the LSP again with nearly a 100% branching fraction. For point C, although 
ti — y bxt is the dominant decay mode, the stop also has sizable branching into tXi, tx\ and 
bxj ■ As for the sbottoms, they have sizable branching fractions into both the top-chargino 
and the bottom-neutralino modes. The former, though slightly suppressed on account of 
phase space considerations, is particularly interesting in that it leads to tops in final states. 
With the stop and sbottom being so light in this scenario, it is obvious that the gluino 
decay branching fractions into stop and sbottom (accompanied by a top or a bottom, as 
the case may be) are significantly enhanced as compared to the typical mSUGRA case. In 
fact, these modes, all of comparable magnitudes, together turn out to be overwhelmingly 
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Decay modes 
(squark/gluino) 


mSUGRA-A 


mSUGRA-B 


mSUGRA-C 




26.6 


27.4 


33.0 


g -» xj*6 


19.0 


19.0 


12.0 





22.0 


22.0 


31.8 




33.0 


35.4 


54.4 


x° -> xf wrf 


32.4 


32.0 


12.0 




10.0 


15.4 


6.0 


X+ -> Xiud 


66.6 


66.6 


66.6 


xt ->■ Xi^i 


33.0 


33.0 


33.0 



Table 6: The branching ratios(%)of the dominant decay modes of the gluinos, lighter neu- 
tralino and chargino states for mSUGRA for the different benchmark points. 

dominant. This, obviously, results in an enhanced scope of having top and/or bottom quarks 
with a high multiplicity. This, in turn, makes it likely to have several leptons in the final 
state (typically from the top quark decays). This particular character of the spectrum, thus, 
raises hopes for A£ (with i = e, /i) final states as a viable signal of SUSY. We, nonetheless, 
do not limit ourselves to these alone, but consider each of the following final states: 

• Opposite sign dilepton (OSD) : {IV) + (> 2) jets + $ T , 

• Same sign dilepton (SSD) : (£ ± f ±) + (> 2) jets + $ T , 

• Trilepton {3£ + jets): 3£ + (> 2) jets + $ T , 

• Hadronically quiet trilepton 4 (3£): 3£ + 1£t , 

• Inclusive 4-lepton (4£ + X): U + X + $ T , 

where I stands for final state electrons and/or muons, denotes missing transverse energy 
and X denotes any associated jet(s). 

Of the various final states listed above, only the hadronically quiet trileptons have their 
origin in electroweak processes such as x^Xi production. However, as can be seen from our 
event selection criteria set down in section 3.2, strong processes which do not give rise to 

4 Thcsc get contributions from electroweak production of a chargino and a neutralino. 
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hard enough jets can also be responsible for such final states. The large rate of tt and bb 
production in NUSM thus leads to relatively higher rates for hadronically quiet trileptons. 
On the whole, rates are never found to exceed a few percent of those with accompanying 
hard jets. 

As is well known, in the LHC environment, even if the hard scattering process were to lead 
to a purely non-hadronic final state, the actual observable final state would, nonetheless, still 
include typically a few jets, originating from underlying events, pile up effects and ISR/FSR. 
In view of this, we define a hadronically quiet event to be one devoid entirely of any jet with 
ET jet > 100 GeV. This avoids unnecessary removal of events accompanied by relatively 
soft jets. 

3.2 Detection and Kinematical Requirements 

Before we mention the selection cuts, we would like to discuss the resolutions of the detectors, 
specifically those applicable to the ECAL, the HCAL and the muon chamber that have been 
incorporated in our analysis [55]. This is particularly important for reconstructing missing- 
Eti which is a key variable for discovering physics beyond the Standard Model. 

We assume that all charged particles with px > 0.5 GeV are detected 5 as long as they 
emanate within the pseudorapidity range \r]\ < 5. For muons though, the applicable pseudo- 
rapidity range is determined by the geometry of the muon chamber to be 6 \rj\ < 2.5. All the 
particles thus detected constitute the "physics objects" that are reconstructed in a collider 
experiments, and are further classified as 

• isolated leptons; 

• hadronic jets formed after identifying isolated leptons; 

• unclustered energy comprised of calorimetric clusters with px > 0.5 GeV (ATLAS) 
and \r]\ < 5, that are not associated with any of the above types of high-Er objects 

Electrons and muons with px > 10 GeV and |?7| < 2.5 may be identified as isolated 
leptons. In order to be deemed isolated, the lepton should be sufficiently separated from 

5 This threshold is specific to ATLAS, while for CMS, px > 1 GeV is applicable. Our results, though, are 

largely insensitive to the exact figure. 

6 Although it seems that muons in the range 2.5 < \rj < 5 would leave their footprints in the tracker, we 

deliberately choose to be consistent with the above criteria. Once again, the inclusion of such muons would 

make little quantitative difference. 
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any other lepton in that it must satisfy ARu > 0.2, where AR = (Arj) 2 + (A0) 2 is the 
separation in the pseudorapidity-azimuthal angle plane. Similarly, it must be far away 
(ARij > 0.4) from all putative jets with Et > 20 GeV. And, finally, the total energy deposit 
from all hadronic activity within a cone of AR < 0.2 around the the lepton axis should be 
< 10 GeV. 

Jets are formed with all the final state particles after removing the isolated leptons from 
the list with PYCELL, the inbuilt cluster routine in Pythia. The detector is assumed to stretch 
over the pseudorapidity range \rj\ < 5 and is segmented into 100 equal-sized (in 77-spread) 
strips. Similarly, the entire 2ir azimuthal spread is again segmented into 64 equal-sized 
strips resulting in a 100 x 64 grid of cells. To register a signal, a minimum Et of 0.5 GeV 
needs to be deposited in a cell, while the minimum Et for a cell to act as a jet initiator is 
assumed to be 2 GeV. All objects within a cone of Ai?=0.4 around the jet initiator cell are 
considered for jet formation, and for a conglomeration to be considered a jet, it must satisfy 
Eobjccts Et > 20 GeV. 

Now, as has been mentioned earlier, all the other final state particles, which are not 
isolated leptons and are yet separated from jets by AR > 0.4 are considered as unclustered 
objects. This includes all electromagnetic objects (muons) with 0.5 GeV < Et < 10 GeV 
and \r)\ < 5(2.5) as well as hadronic energy deposits with 0.5 GeV < Et < 20 GeV and 
1 7/| < 5. Such unclustered energy deposits need to be taken into account in order to properly 
reconstruct any missing- Et- 

Any detector suffers from finite resolutions and collider detectors are no exception. To 
approximate the attendant experimental effects, we smear the energies (transverse momenta) 
with Gaussian functions. Nominally, the widths of the distributions have different contri- 
butions (accruing from different sources), each with a characteristic energy dependence and 
with these being added in quadrature. To wit (all energies are measured in units of GeV), 

• electron/photon energy resolution: 

u(E) cl c 



E J~E E 



where 



(6a) 



(0.030,0.005,0.2) H<1.5 
(a, 0, c) = < (ob) 



(0.055,0.005,0.6) 1.5 < ] ^7 ] < 5 



muon pt resolution : 




a(P T ) J a \rj\< 1.5 

Mog ToHeV l-«<N<2-5 
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with 



(a,b) 



(0.008,0.037) 
(0.020,0.050) 



\r)\ < 1.5 

1.5 < \t]\ < 2.5 



(7b) 



• jet energy resolution : 



<t(Et) a 



(8) 



with a = 0.55 being the default value used in PYCELL 



• unclustered energy resolution : 



a 



(E T ) = a^E? nc - E ^ 



(9) 



where a ~ 0.55. One should keep in mind here that the x- and y-components of 
j^Unc.E neec [ sm eared independently (with identical widths). 

Once we have identified the 'physics objects' as described above, we sum vectorially the 
transverse components of all the momenta smeared thus to obtain the total visible transverse 
momentum. Clearly, the missing transverse energy is nothing but the magnitude of the visible 
transverse momentum, namely 



where the sum goes over all the isolated leptons, the jets as well as the unclustered energy 
deposits. At this stage, we are in a position to impose the selection cuts, namely 

• Missing transverse energy Jf}? > 100 GeV, 

• p l T > 20 GeV for all isolated leptons, 

• E^ > 100 GeV and \r) jet \ < 2.5, 

• For the hadronically quiet trilepton events, as also for inclusive A£ events, we reject, 
in addition, any event with a same flavour opposite sign lepton pair satisfying \Mz — 
Mi+£- 1 < 10 GeV. Such events are characterised by the requirement of having no 
central jet with E T > 100 GeV. 

We have generated the corresponding SM backgrounds (with identical kinematical cuts) 
using Pythia. The bulk of the contribution comes from ti events. To take into account 




(10) 
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the next to leading order (NLO) and next to leading log resummed (NLL) corrections — not 
included in Pythia — we rescale the results by the appropriate i^-factor [56] viz. 2.23. Ex- 
clusive diboson (WW, WZ, ZZ) production constitutes another potential background, but 
it is easy to see that except for the hadronically quite trilepton channel, these contributions 
are very sub- dominant. Furthermore, these are reduced drastically by the cuts imposed, 
especially by the one on the leptonic invariant mass. Inclusive, i.e. including (multi-)jets, 
gauge boson production is another very serious background, but can be estimated with a 
high accuracy using Alpgen [57]. The combination of a large missing along with the 
requirement of at least two high-p^ leptons reduces even this to innocuous levels. 

3.3 Results 

The event rates in the various channels discussed in the preceding section would, of course, 
differ amongst themselves and also depend on the point in the parameter space, both on 
account of the differing production cross sections and branching fractions as well as the 
kinematical restrictions imposed. In Table 7, we tabulate the event rates in different chan- 
nels obtained for the points A, B and C of the NUSM scenario as well as those for the 
corresponding mSUGRA ones. Also shown are the respective SM backgrounds. 

For the NUSM benchmark points, the gluino decays dominantly into top-stop and bottom- 
sbottom (see table 4). The source of leptons in the final state can thus be both the stop and 
the sbottom which can lead to the top and chargino in the next stage of the cascade. Of 
course, appropriate branching and combinatoric factors are to be used in each case. For the 
mSUGRA benchmark points, on the other hand, the gluino (which is the lightest strongly 
interacting superparticle) decays primarily into the three-body channels such as ttXi and 
tbxf (see table 5). It should be remembered, however, that gluino decays mediated by light 
squark flavours are not entirely negligible and in fact, can account for upto one-third of 
the decays. This causes an effective enhancement, in the NUSM cases, of decays into the 
intermediate states containg top/chargino. 

Another crucial difference is the splitting between the Xi an d Xi states. With the splitting 
being large in the NUSM case, the W ± from Xi can be nearly on-shell (as opposed to an 
off-shell one in the mSUGRA case), thereby resulting in typically harder leptons. Together, 
these features are responsible for effectively reducing the rates for leptonic final states for 
mSUGRA in comparison with the NUSM benchmark points consistent with the dark matter 
constraints. Of course, the already mentioned difference in gluino masses has also a small 
role to play. 
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Model Points 


(JOSD 


&SSD 






0"4£ 


A 


103 


24.0 


14.9 


3.1 


3.1 


mSUGRA-A 


33.7 


15.4 


8.1 


0.4 


1.3 


B 


135 


28.7 


19.0 


4.4 


3.8 


mSUGRA-B 


38.9 


16.9 


9.1 


0.4 


1.5 


C 


23.9 


7.3 


2.9 


0.1 


0.3 


mSUGRA-C 


1.8 


0.6 


0.4 


0.1 


0.1 


SM Backgrounds 


tt 


1.10 x 10 3 


18.1 


2.7 


5.3 


0.0 


ZZ,WZ,ZH,Z~i 


16.3 


0.3 


0.5 


1.1 


0.4 


Total Bkgd 


1.12 x 10 3 


18.4 


3.2 


6.4 


0.4 



Table 7: Event-rates (fb ) after cuts for non-universal scalar mass points and corresponding 
mSUGRA points with same gluino mass. CTEQ5L parton distribution functions used with 
[If = Vr= V§- 

Based on the above observations, the following features in the results are noted. 

• For each of OSD, 3£ and 4£ final states, the difference in the absolute rates between 
the NUSM points and the corresponding mSUGRA ones is remarkably large. As has 
been argued at the beginning of Section 3.1, this can be understood in terms of the 
relative lightness of the third-generation squarks in the NUSM scenario. While it might 
be tempting to aver that this alone would serve to distinguish NUSM from mSUGRA 
scenarios, a little reflection shows that just the absolute rates are not enough for this 
purpose and a combination of observables would be required. 

• For the same sign dilepton final state (a manifestation of the Majorana nature of the 
gauginos), the signal to background ratio (S/B) exceeds unity for both points A and 
B. For point C, while S/B < 1, a discovery is possible with an integrated luminosity 
of only 10 fb -1 . This, though, is not surprising, for SSD is well known for its efficacy 
in SUSY search. Note though that the SSD rate cannot really distinguish between 
NUSM and its mSUGRA counterpart. While the case for point C may look promising, 
it is precisely the case where the rate (and, hence, statistical significance) is low. 

• While S/B^l for the 4-lepton final state as well, the smaller rates for this signal sig- 
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nificantly reduces its potential as a discovery channel. However, this could potentially 
serve as a very efficient discriminator between scenarios. For example, the S/B ratio 
is close to unity even for point C. 

• The situation for the 3£ + jets is somewhat better than the 4£ one. The rates are larger 
while maintaining the difference between NUSM and mSUGRA. Once again, for point 
C, the S/B ratio is close to unity, while for mSUGRA-C, it is suppressed. 

• As for the OSD and the hadronically-quiet 3£ final states, generically, S/B < 1. The 
former though boasts of the largest event rates. For points A and B, this signal begins 
to stand well over the background fluctuation for an integrated luminosity of as little 
as 1 fb" 1 whereas 2 fb" 1 would allow a discovery claim. For the hadronically-quiet 
trilepton mode, the required luminosity is ~ 10fb _1 for points A and B. Qualitatively, 
these two points are very similar to each other, especially as far as the superpartner 
masses are concerned. The main difference lies in the Higgs sector, which has not been 
explored here. 

• For point C (and even more so for mSUGRA-C), the heaviness of the spectrum trans- 
lates to lower event rates. 

• It is obvious that, for an integrated luminosity of 30 fb _1 , many of these channels would 
allow us to register a 5a discovery claim 7 . In Table 8, we summarise this information 
for each of the channels and parameter points, both in the NUSM scenario as well as 
their mSUGRA counterparts. 

We now discuss the profile of the dominant (though not necessarily the most background- 
free) signal mode, namely events with opposite sign dileptons. The quest is to see if quanti- 
tative features in the same could be used to either accentuate the discovery potential or as 
discriminators between models and/or parameter points. In Fig.2, we display the normalized 
(to unity) distributions of missing transverse energy, a most crucial aspect of supersymmetry 
signals. As a comparison of the first two panels shows, the Ifa distributions for parameter 
points A and B look very similar, which is but a consequence of the aforementioned simi- 
larity between the corresponding spectra. Furthermore, all of them are discernibly different 
from those for the corresponding mSUGRA points 8 . That the latter are softer can be under- 

7 The required luminosity is much smaller for some channels and parameter points. 

8 Note that we are concerned here about the shape of the curve, not the absolute magnitude, which, of 
course, are different (vide Table 7.) 
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mSUGRA-A 
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V 


B 


V 


V 


V 


V 


V 


mSUGRA-B 


V 


V 


V 


X 


V 


C 


V 


V 


V 


X 


X 


mSUGRA-C 


X 


X 


X 


X 


X 



Table 8: ha visibility of various signals for an integrated luminosity o/30fb . A y/ indicates 
a positive conclusion while a x indicates a negative one. 

stood by realizing that the main production channel for the mSUGRA spectrum is pp — >■ gg 
and that, unlike in the NUSM case, the gluino undergoes a three body decay, resulting 
in relatively less momentum imparted to the LSP. Note also that the the dominant (tt) 
background — as displayed in the first panel — is almost as soft as the mSUGRA signals, and 
thus a hardening of the IpT cut would have considerably improved the S/B ratio for the 
NUSM cases, while worsening it for the mSUGRA ones. And, finally, for point C, the dif- 
ference is even more stark. However, with the size of the signal being small in this case, it 
is not immediately apparent whether this could be used to any advantage. 

Another kinematical variable often used advantageously in searches for new physics in- 
volving JpT is the "effective mass" defined to be the scalar sum of the transverse momenta 
of the isolated leptons and jets and the missing transverse energy, viz. 

m cS = ^(p^iso.i + ^(P^jets + $T ■ (11) 

In Fig.3, we display the corresponding distributions, again for both the NUSM points and 
their mSUGRA counterparts. As the first panel shows, as far as point A is concerned, there 
is little to choose between this distribution and that for the corresponding mSUGRA point. 
Similar is the case for point B (second panel). Note, furthermore, that the peak in either 
case is at a fairly large value of m e ff. While this, at first sight, might seem contradictory to 
the oft-repeated claim that this distribution should peak roughly at twice the mass of the 
dominant particle being produced, the reason for this discrepancy is easy to appreciate. First 
and foremost, with the strong demands made on the transverse momenta of the two leading 
jets, the contribution from stop-pair production reduces drastically. This is understandable 
since the relatively small difference between the stop mass and those of the lighter chargino 
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Missing Et (GeV) Missing Et (GeV) " Missing Et (GeV) 

Figure 2: Missing transverse energy distribution (normalized to unit area) for opposite sign 
dilepton (OSD) events. The eponymous panels refer to the respective representative points 
in the parameter space. Also shown are the analogous distributions for the corresponding 
mSUGRA points. The first panel also shows the distribution accruing from the overwhelm- 
ingly dominant SM background, namely ti production. 

implies that the b from stop decay tends to be softer. With the stop-pair contribution 
thus being effectively decimated 9 , this also offers hints as to why the NUSM and mSUGRA 
distributions look so similar. And, with the gluinos themselves being produced with a 
considerable transverse momentum, it is easy to understand why the distribution peaks at 
a high value of m e s. 

For point C though, a remarkable difference with the mSUGRA counterpart is immedi- 
ately apparent. This can be traced to the fact that the charginos are considerably softer for 
mSUGRA-C as compared to the point C. This allows for the second, subdominant, hump 
slightly above twice the mass of the chargino. While this could, in principle, be used to dis- 
criminate between the two scenarios, unfortunately the smaller rates tend to make the task 
a difficult one. And, finally, note that the distribution for the (dominant) ti background is 
considerably softer than those for any of the six supersymmetric parameter points discussed 

9 It might seem paradoxical that we are altogether sacrificing the signal from the light stop, a cornerstone 
of this scenario. However, including the stop contribution would require softening the pt requirements, a 
process fraught with danger in the context of the LHC. In the absence of a full-scale simulation including 
multiple scattering and underlying events, we deliberately desist from this. 
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Figure 3: As in Fig. 2, but for the effective mass instead. 



here. Thus, imposing a requirement such as m c g ^ 2 TeV would have significantly improved 
the S/B ratio for the OSD signal. This, though, would have eliminated the secondary hump 
for the mSUGRA-C case. 

Having discussed the prospects of refining and/or using the kinematical variables in the 
OSD sample towards discriminating between scenarios, we now consider a set of observables, 
namely the ratios of events seen in various channels. As is well known, there is a great 
advantage to the use of such variables in that it almost entirely eliminates some systematic 
uncertainties such as that in the luminosity and drastically reduces others such as those 
corresponding to the choice of the parton distributions, the choice of renormalization and 
fragmentation scales etc. In Fig. 4, we present the ratio of the subordinate channels with the 
dominant (OSD) channel for each of the parameter points. 

At this stage, we can easily formulate the means of discriminating between a NUSM 
point and the corresponding mSUGRA one, namely 

For parameter points (such as A and B) with a relatively smaller m 1( / 2 but large m 
(i.e., when the gluino is considerably heavier than the stop/sbottom but sufficiently 
lighter than the first two generation squarks and sleptons), the NUSM scenario would 
typically result in a smaller proportion of same sign dilepton events, as is clear from 
FigA(a), when compared to the OSD rates. This can be attributed to the fact that 
the OS'.D-rate increases significantly for the gluinos decaying through stop-top (with, 
consequently, OSD being possible even from the decay of a single gluino, irrespective 



25 




Benchmark Points Benchmark Points 



Figure 4: Event ratios with respect to opposite-sign dilepton(OSD) events for NUSM and 
mSUGRA cases at the points A,B,C mentioned in the text. Colour Code: Red: mSUGRA, 
Green: NUSM scenario. 

of how the other one decays), whereas the SSD relies on the good old fact of the 
gluino being a Majorana spinor, with only a slight increment to the leptonic branching 
fraction due to the decays through third generation. 

"A" The mSUGRA-A and mSUGRA-B points have sufficiently small values of \i and these 
points indeed fall in the HB/FP region. This implies that there is more of Higgsino 
in the lighter chargino (xf ) and the second lightest neutralino (x^)- Therefore, the 
leptonic signals are weakened compared to what we predict in NUSM. 

"A" Again, for points with a small m^j but large m (such as A and B), the rates (absolute 
and relative) for the hadronically quiet trilepton mode are markedly higher for the 
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NUSM case. This can be attributed to the aforementioned feature of the NUSM 
spectrum which renders it easier to have large-pr isolated leptons. The turnaround 
for point C vs mSUGRA-C is a consequence of the lightness of the charginos in the 
latter allowing a large contribution through the production of the charginos and the 
second-lightest neutralino. 

~k Overall, it is self-evident that a combination of these ratios would serve to easily 
distinguish between the two scenarios. 

it And, finally, the relative lightness of the stop and the sbottom (and the consequent 
fact of the gluinos decaying through these), renders the NUSM signal b- rich. Invoking 
6-tagging (which we had not done in the results presented so far) would thus present 
us with a very useful discriminator. With this in view, we perform a study in OSD 
channel associated with two or more partonic fe-jets (>2b + OSD). We assume a b- 
tagging efficiency [58] of e& = 0.5 for pt > 40 GeV and \r)\ < 2.5. The OSD event 
selection criteria remain the same. As expected, we see a clear distinction between the 
NUSM and the corresponding mSUGRA ones in the absolute event rates as shown in 
Table 9. With the NUSM sample being particularly rich in 6's, the suppression in rates 
as compared to those in Fig. 7 is understandably less severe than ef. This, of course, 
does not apply as well to point C and far less to the MSSM cases, resulting in a large 
suppression of the latter. 



Model Points 


A 


mSUGRA-A 


B 


mSUGRA-B 


C 


mSUGRA-C 


tt 


>2b + OSD 


36.6 


6.3 


46.4 


10.2 


6.7 


0.5 


148.7 



Table 9: Event rates (fb) at different benchmark points and for the tt background for a final 
state >2b + OSD. 

4 Summary and Conclusions 

We have studied a case of nonuniversal scalar masses, wherein the first two families of 
squarks as well as sleptons of all generations are much heavier than the third family of 
squarks and the Higgs scalars. The universality of gaugino masses has been adhered to. 
We confine ourselves only to that region of the parameter space where one achieves a relic 
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density consistent with the WMAP data. LSP annihilation is efficiently mediated by the 
pseudoscalar Higgs, with the 'funnel region' being significantly extended toward small values 
of tan/3 when compared to mSUGRA. Having ensured that the region of SUSY parameter 
space thus isolated is consistent with all constraints from FCNC and CP-violation, we have 
proceeded to investigate the signals of this scenario at the LHC. 

Although stop-pair production is the dominant SUSY process in such scenarios, the 
stringent cuts that we choose to impose results in stop-cascades being suppressed. Rather, 
the dominant contribution to the signal rates turns out to be gluino pair-production. The 
relatively large multiplicity of top quarks produced in the cascades results in enhanced rates 
for two, three and four-lepton final states, together with missing-E^ and hard jets. In 
particular, the usefulness of four-lepton final states is highlighted through this analysis. 

Based on the study of a few benchmark points (corresponding to different gluino masses 
and tan j3), we find that, using the criteria chosen by us, it is possible to probe the above sce- 
nario with an integrated luminosity of 30 fb" 1 for gluino masses up to about 1.2 TeV. (Indeed, 
for certain significant parts of the parameter space, even 2 fb _1 would be enough.) The reach 
can be potentially extended further once more luminosity accrues. We also demonstrate that 
it is possible to distinguish this scenario from an mSUGRA-one tuned at the same gluino 
mass and satisfying the WMAP constraints. The usefulness of the ratios of events in various 
channels is clearly elicited from our study. Moreover, such distinction is facilitated by the 
effective mass distribution of events for gluino masses on the higher side, i.e. above a TeV. 
Thus, we succeed in illustrating that a multichannel analysis is not only able to probe such 
nonuniversal SUSY scenarios satisfying the relic density constraints, but can also highlight 
notable differences with a simple-minded model based on universal SUGRA. 
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